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Continuum Widom-Rowlinson model: In the beginning of seventies, a great 
progress has already been made in understanding equilibrium phase transitions 
and critical phenomena in classical (i.e. non-quantum) lattice models, where ba- 
sic elements are localized on vertices of a regular discrete structure and interact 
with their nearest neighbors. The lattice gas analogy [1] of Onsager's exact solu- 
tion of the two-dimensional Ising model enabled one to understand qualitatively 
as well as quantitatively physical implications of the spontaneous breaking of the 
particle-hole symmetry on the existence and critical properties of high-density liq- 
uid and low-density vapor phases. On the other hand, little was known about the 
liquid- vapor equilibrium in continuum models (fluids), where the distance between 
interacting particles is a continuous variable. In 1970 Widom and Rowlinson (WR) 
introduced a simple model of the classical fluid in thermodynamic equilibrium [2], 
found a counterpart of the particle-hole symmetry and described implications of 
this symmetry on a potential liquid-vapor phase diagram. Subsequent studies 
proved rigorously the existence of liquid- vapor phase transitions in the WR and 
related models in spatial dimensions d > 2; up to now WR fluids are the only con- 
tinuum systems with decaying interactions for which on has this kind of rigorous 
results. 

The WR model consists of identical molecules living in an infinite <i-dimensional 
continuous space of points r 6 V (V — > R d ), around the center of each molecule 
there is a microscopic sphere of radius a and volume vq. The potential energy U 
associated with a given configuration r 1 , . . . ,r N of N molecules is defined by 

Ufa, ...,r N ) = [Wfa, ...,r N )- Nv ] e/v , (1) 

where W denotes the volume covered by the corresponding N (in general interpen- 
etrating) spheres and e > is some energy constant. Since W fulfills the evident 
inequalities v < Wfa, . . . , r N ) < Nv , the potential energy is bounded as follows 

-{N-l)e<Ufa,...,r N )<0; (2) 

the lower bound ensures a correct thermodynamics, the upper bound tells us that 
the short-ranged forces among molecules are purely attractive. The model is usu- 
ally studied in the grand canonical formalism characterized by the dimensionless 
inverse temperature 9 = e/ (k B T) and particle activity z, the corresponding particle 
density (N)/V and pressure P are considered in the dimensionless combinations 

p(z,6)=v (N)/V, V(z,6)=Pvo/(k B T); (3) 
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z is normalized so as to be asymptotically equal to p in the ideal gas limit p — > 0. 

The symmetry of the WR fluid, whose spontaneous breaking is responsible for 
the existence of the liquid and vapor phases, is hidden in the original formulation. 
It becomes transparent after mapping the WR model in a thermodynamic sense 
onto a two-component WR mixture with non-additive hard-core interactions. 

The WR mixture is composed of two kinds of molecules a G {A, B} interacting 
pairwisely U({r}) = Y^i< 3 <v>ap(\ti ~ v j\)i where the particles of the same species 
do not interact whereas the unlike species interact with a hard-core repulsion at 
distances smaller than 2a, i.e. 

_ f oo if a ^ (3 and r < 2a, , . 

Ua/3[T) -\ otherwise. 1 > 

Let za and z B be the dimensionless activities of the two species, normalized so as to 
be equal to the densities pA = v (Na) /V and p B = v (N B ) /V in the ideal gas limit, 
and V m ix(zA, z B ) be the dimensionless pressure (independent of the temperature 
since the only interactions are infinitely strong repulsions) defined in analogy with 
(J3]). The mixture possesses the obvious symmetry with respect to the interchange 
of A and B molecules, P m ix(^A, z B ) = "Pmix^B, za)- Let us consider the line of 
equivalent activities za = zb = z. The system behaves like the ideal gas at very 
low z, i.e. there is just one pure (mixed) phase with equivalent particle densities 
Pa = Pb- At very large z, since unlike molecules experience an infinitely strong 
repulsion, a mixed phase suffers from packing effects, which are substantionally 
reduced in a demixed phase with a single majority component. Consequently, the 
A — B symmetry of the WR mixture is broken and the system can exist in two 
different homogeneous (i.e. translation invariant) pure phases: the A-rich phase 
with pa — Pb = 8p > and the complementary S-rich phase with ps — Pa = 8 p. 
The two phases become equivalent (Sp = 0) at the "demixing" critical point z^. 

The two-component WR mixture is exactly solvable in one dimension with no 
phase transition. For dimensions d > 2, the proof of the existence of more than 
one pure thermodynamic phase for sufficiently large z was given by Ruelle [3j using 
the Peierls contour method. 

Theorem 1. For each dimension d > 2, there exists Zd > such that the A-B 
WR model on the line of equivalent activities za = zb = z has 

• A unique homogeneous mixed phase for all z < Zd- 

• Multiple homogeneous demixed phases for all z > Zd- 

By integrating over the coordinates of one species in the two-component WR 
mixture, its thermodynamic equivalence with the original WR model is realized 
via the identity between the dimensionless pressures [2j 

V(z,6) = -e + V aSx (6,ze°). (5) 

The line za = zb of the A-B mixture is thus transcribed to the line of symmetry 

z = 8exp(-6) (6) 
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in the original WR model and the demixing critical point Zd has its image at the 
critical point (9* = Zd,z* = Zde~ Zd ). For 9 > 9*, the line of symmetry ([6]) is the 
coexistence curve of the gas (vapor) phase with particle density p g and of the liquid 
phase with density pi > p g . The mapping enables one to deduce, without solving 
explicitly the WR model, singular behaviors of thermodynamic quantities on the 
line of symmetry ([6]) in the neighborhood of the demixing critical point. Some 
features of singular behaviors differ from those predicted by the lattice gas model 
PQ, e.g. the temperature derivative of the mean density (p g + pi)/2 is proportional 
to the heat capacity at constant volume and goes to — oo at the critical point [2]. 

The multicomponent generalization of the A-B WR mixture consists in consid- 
ering molecules of M different types a = 1, 2, ... , M, all having the same activity 
z. The molecules interact pairwisely, the only interaction is a hard-sphere repulsion 
between any two particles of unlike species (JH). 

It was shown [I] that for any finite number M of components the WR model in 
d > 2 dimensions exhibits the demixing phase transition at some Zd{M) (Theorem 
1); in a pure demixed phase, the homogeneous density of just one of the components 
is dominant, say p\ > p 2 = P3 = . . . = pu- 

A version of the M-component WR model, in which the hard-sphere interaction 
between unlike species is slightly modified to that of parallel hard (hyper-)cubes, 
was studied in the limit of infinite dimensionality d — > oo [5]. In that limit, the 
calculation of thermodynamic functions within the second virial coefficient is exact. 
It turns out that for M > 31 the transition from the mixed phase at small values 
of z to the demixed phase at large values of z is preempted by solidification at 
intermediate values of z, z c (M) < z < Zd(M). In the corresponding crystal phase, 
all species are equivalent (pi — p 2 — . . . — pu — p/M with p being the total 
density of molecules), but the density p = p(r) varies periodically in space, i.e. the 
translation symmetry is broken. The origin of this phenomenon is purely entropic: 
for a large number of component M it pays the system to create a periodic structure 
of alternating dense and sparse regions, where the particles in the dense regions 
are less restricted by the hard-core repulsions coming from particles in the sparse 
regions. 

Lattice Widom-Rowlinson model: The lattice version of the two-species WR 
model was introduced by Lebowitz and Gallavotti [6j. The model is defined on 
a regular ^-dimensional integer lattice Z d , each lattice site % can be either empty 
{local state cxj = 0, the corresponding activity z(0) = 1} or singly occupied by a 
particle of type A {<7j = 1, z(l) = z} or B {a, = 2, z(2) = z}. The only interaction 
is an infinite repulsion between unlike particles on nearest-neighbor pairs of lattice 
sites, i.e. the potential energy associated with a given state configuration a of 
all sites is formally expressible as U(cr) = Y,(i,j) u {. (y ii°'j)i where the interaction 
potential between nearest-neighbor pair of lattice sites is given by 

J oo if <n ^ Oj and a { £ 0, a s £ 0, , , 

^^'J-jo otherwise. UJ 

A relatively simple proof of the demixing phase transition (Theorem 1) in the lattice 
WR model was done using standard Peierls methods [6]. Moreover, in dimensions 
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d > 3, the existence of a sharp interface between coexisting demixed A-rich and 
B-rich phases, generated via the respective fixed A and B boundary conditions on 
opposite sides of a box, was proved in the limit of the box size going to infinity [7j. 

In the multicomponent lattice generalization of the two- component WR model, 
each lattice site % G Z d can be either empty {crj = 0, z(0) = 1} or singly occupied 
by a particle of type <Ji = 1, 2, . . . , M, all having the same activity z. The potential 
energy of a state configuration a reads U(a) = u(a^ crj), where the interaction 
potential between nearest-neighbor sites is given by relation (EJ). Note that 
if one replaces oo in (J7J) by some finite U ^ 0, the system is equivalent to a dilute 
Potts model. The number density of species a — 1, . . . , M at site i will be denoted 
by Pi(cr), the total density of particles at site i by pi — J2cr=i Pi{ a )- 

The multicomponent lattice WR model was studied by Runnels and Lebowitz 
[I]. In dimensions d > 2 and for any finite number M of components, the lattice 
WR model exhibits the demixing phase transition at some z d (M) (Theorem 1); in 
a pure demixed phase, the site-independent density of just one of the components 
is dominant, say p(l) > p(2) = . . . = p(M). When the number of components M 
is larger than some minimum M , an entropy-driven crystal phase prevents from a 
direct transition between the mixed and demixed phases. 

Theorem 2. For each dimension d > 2, if the number M of components of the 
lattice WR model exceeds some M , there exist two activity thresholds, the crystal 
one z c (M) > and the demixed one z d (M) > z c (M), such that: 

• For z < z c (M), there exists a unique translation invariant mixed phase with 
species- independent densities Pi(cr) = p/M. 

• For intermediate values of z, z c {M) < z < Zd(M), there exist two distinct 
crystalline phases in which one of the alternating sublattices (even or odd) 
is preferentially occupied by particles, i.e. the translation symmetry (but 
not the species symmetry) is broken. More precisely, the average particle 
densities on the even and odd sublattices are unequal, p e ^ p Q , while the 
average densities of the species a = 1, 2, . . . , M are the same within a given 
sublattice, i.e. p e (cr) = p e /M an p (a) = p /M. 

• For z > Zd{M), there exist M translation invariant demixed phases with the 
density dominance of just one of the species, say p(l) > p(2) = . . . = p(M). 

The transition at z c (M) is of second order. For asymptotically large M the crystal 
threshold z c (M) behaves like A cr /M, where A cr is the critical activity of the one- 
component lattice gas with nearest-neighbor hard-core exclusions. The rigorous 
upper bound M < 27 6 derived for the square lattice [1] was surprisingly large and 
made the crystallization phenomenon only of academic interest. 

The exact solution of the multicomponent WR model on the Bethe lattice of 
coordination number q gives M = \qj (q — 2)] 2 [8j, which would suggest more real- 
istic Mo ~ 4 for the square lattice with q = 4. The Monte-Carlo (MC) simulations 
for the square lattice [8] imply M = 7. The MC simulations indicate that the tran- 
sition at Zd{M) is of second order for M < 4 and of first order for M > 5, putting 
the M-component lattice WR model in the M-state Potts model universality class. 
For asymptotically large M, z d (M) ~M-2 + l/MH . 
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The consideration of hard-core exclusion to next-to-nearest-neighbors leads to 
analogous phases which numbers and characters depend on the specific geometry 
of the model [9j. 
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